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Abstract. Statistical arguments show that the volume- and time-averaged kinetic 
power of GRBs and fireball transients (FTs) into an L* galaxy like the Milky Way 
is at the level of 10 40 ergs s^ 1 . This number, though with wide uncertainties related 
to the internal or external shock efficiency, is sufficient to power hadronic cosmic rays 
observed locally. 

The release of energy by the high-mass progenitor stars of GRBs and FTs is sufficient 
to power the ultrahigh energy cosmic rays, as already shown by Waxman and Vietri 
in 1995. It is sufficient to power the cosmic rays above the knee of the cosmic-ray 
spectrum. Indeed, all hadronic cosmic rays could originate from the high-mass ( ^ 100 
Mq) stars that collapse to black holes, in the process forming GRBs and FTs. This 
source class represents a new solution to the problem of cosmic-ray origin. 

The ~ 10 4 -10 7 black holes made by these stars could make their presence known by 
radiating as they accrete from the ISM, by microlensing background radiations, and by 
forming luminous binary systems. Some unidentified EGRET sources could be isolated 
black holes that accrete from the ISM. Better imaging and sensitivity with GLAST 
and TeV observatories will test this model for the unidentified 7-ray sources, and this 
theory for cosmic-ray origin. 



I INTRODUCTION 

Fireball transients are those explosive events that propel a significant fraction of 
their ejecta kinetic energy in the form of relativistic baryonic outflows. FTs include 
dirty (r ^ 300), loaded (r ~ 300), and clean fireball (r ~ 300) subclasses [1], of 
which BATSE is most sensitive to those explosions with initial bulk Lorentz factors 
T ~ 300. The dirty and clean fireball classes of exploding stars are difficult to 
discover because of design and sensitivity limitations of detectors flown to date [1,2], 
though the X-ray flashes detected with Beppo-SAX [4] could be representatives of 
the dirty fireball class. HETE-II, Swift, and GLAST can all be expected to make 
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progress on this front, though a lobster-eye X-ray telescope seems most promising 
for discovering dirty fireballs. 

In a recent paper [3], I deduced the time- and volume-averaged kinetic power of 
GRBs and fireball transients (FTs) in the Milky Way on the basis of a statistical 
treatment employing the external shock model of GRBs [5,6] for both the prompt 
7-ray luminous and afterglow phase, performed by Bottcher and myself [7]. There 
I found that the FT power into the galaxy was within an order of magnitude of 
the inferred average Galactic cosmic-ray power, and I considered the possibility 
that the high-mass stellar progenitors of FTs were the sources of the cosmic rays 
(CRs) above and below the knee of the CR spectrum. I also discussed observational 
difficulties in the scenario advanced by Ginzburg and Syrovatskii [8] that hadronic 
CRs are powered when supernovae collapse to neutron stars. 

Although 7-ray astronomy was predicted to solve the CR origin problem (see [9] 
for a recent assessment), evidence for CR source origin still rests primarily upon 
arguments about available sources of power, which in turn hinge upon statistical 
analyses. The numerical simulation of GRB statistics [7] that was used to derive 
the FT power in the Milky Way modeled detector response and generic GRB spec- 
tral and temporal behavior to fit jointly large samples of GRB properties measured 
with BATSE, namely the peak photon energies e° of the prompt vF v spectra, the 
peak count rates, and the £ 1 s £50 durations. The predicted redshift distribu- 
tion of GRBs can be used to test the model, but will require several score GRB 
redshifts from HETE-II and/or Swift. The redshift sample should preferably come 
from a single mission in order to minimize triggering differences and calibration 
uncertainties between detectors. 

Rather than review the numerical results for a third time (see [10,11]), it seems 
more fruitful to devote this space to an analytic statistical model that confirms in 
all respects the numerical results. We gain clarity at the expense of the accuracy 
found in [7]. 

In section II, the standard argument that only supernova explosions are fre- 
quent and energetic enough to power the galactic cosmic rays is recited, while 
adumbrating certain weaknesses in the underlying assumptions. The statistics of 
cosmological GRBs are treated in Section III. Depending on shock efficiency dur- 
ing the prompt phase, the power of the high-mass progenitors into an L* galaxy is 
evaluated. GRBs and FTs are shown to be capable of injecting sufficient energy 
into the ISM to power the hadronic CRs, as discussed in Section IV. Potentially 
observable consequences of the 10 4 -10 7 black holes that are formed from FTs in 
the Galaxy are considered in Section V, including a possible solution to the origin 
of some of the unidentified EGRET sources. The lineage of stars that explode as 
different novae types is sketched in Section VI. Section VII provides a summary 
and a brief discussion of the acceleration and adiabatic loss problem, that will be 
treated in more detail elsewhere [12]. 



II COSMIC RAY AND SUPERNOVA POWER 



Observations show that TeV electrons are accelerated by SN remnants (SNRs), 
but there is as yet no direct observational evidence for hadronic CR acceleration 
by SNRs [3,14]. Perhaps the most compelling argument that hadronic CRs are 
powered by SNRs is the claim that only SNe inject sufficient power into the Galaxy 
to provide the observed CR energy density [8,13]. The local energy density of 
CRs is ucr ~ 1 eV cm" 3 fa 10~ 12 ergs cm" 3 . The required CR power is thus 
Lcr ~ u CRYgai/tesc, where Vg a i is the effective volume of the Galaxy from which 
CRs escape on a timescale t esc . If CRs are trapped in a disk of 15 kpc radius and 
100 pc scale height, then V ga i ~ 4 x 10 66 cm 3 . Observations of the light elements 
Li, Be, and B that are formed through spallation of C, O, and N indicate that 
CRs with energies of a few GeV per nucleon — which carry the bulk of the CR 
power — pass through fa 10 gm cm" 2 before escaping from the disk of the Galaxy. 
A mean disk density of one H atom cm -3 gives t esc ~ 6 x 10 6 yr, implying that 
Lqr ~ 2 x 10 40 ergs s -1 . Analysis of the composition of isotopic CR 10 Be yields a 
larger value of t esc , implying a smaller mean matter density but a larger effective 
trapping volume of the Galaxy so that, in either case, 

Lcr « « 5 x 10 40 ergs s" 1 . (1) 

tesc 

The galactic SN luminosity L SN ~ (1 SN/30 yrs) x 10 51 ergs/SN w 10 42 ergs s" 1 
which, even given a 10% efficiency for converting the directed kinetic energy of SNe 
into CRs that seems feasible through the shock Fermi mechanism, is completely 
adequate to power the hadronic cosmic radiation. 

Although 7-ray astronomy was supposed to solve the cosmic-ray origin problem, 
this has not happened. The predicted ir° decay feature at 70 MeV has not been 
detected from SNRs by EGRET [15], the Whipple imaging air Cherenkov tele- 
scope has not detected emission consistent with hadronic CR acceleration by SNRs 
[16], and the spectrum of the diffuse galactic 7-ray background contradicts the as- 
sumption that CR protons are uniformly distributed throughout the Galaxy with 
a spectral shape that is the same as observed locally [17]. Moreover, it is becoming 
increasingly clear that the stochastic nature of explosive phenomena in the Galaxy 
is important for interpreting radiation emitted by Galactic CRs [18,19]. 

The EGRET observations of the diffuse galactic 7 radiation voids the assumption 
that hadronic CRs uniformly inhabit the Galaxy, so that the required power could 
be overestimated if CR leptons emit most of the diffuse galactic radiation. Alter- 
nately, we could live in a region of enhanced CR hadron energy density compared 
to the Galactic average. This is more likely if CRs are produced by rare powerful 
events, which then the power requirements are also reduced. Our location in the 
Gould belt shows that we live near a region of enhanced stellar formation activity 
[20]. 

Even if equation (1) is reliable, the massive stars that collapse to black holes 
while making GRBs can provide a kinetic power into the Galaxy comparable to 




that from SNe, as we now demonstrate. 



Ill STATISTICS OF COSMOLOGICAL GRBS 

The new treatment of GRB statistics presented here is motivated by the recent 
Hipparchos, SN la and BOOMERANG results which indicate that we are living in 
a universe with non-zero cosmological constant [21], and by Beppo-SAX and follow- 
on observations demonstrating that GRBs are associated with events occurring in 
star-forming regions. For background on cosmological statistics, see [22-26,7]. 

The rate N(> <p p ) at which observers detect events above some threshold peak 
flux (f) p is given by 

Atic f°° 

N (> <t> P ) = -ft dz ^cosmoW n com {z)P tr [(t> p {z)\ , (2) 

-no ■'0 

where H = 3.24 x 10~ 18 h s _1 , n com (z) is the comoving rate density (cm -3 s _1 ) of 
explosions at redshift z, and <p p (z) relates the observed peak uF u flux to the peak 
source luminosity in the waveband at which the detector is most sensitive. 

The cosmological distance d cosmo = dc(l + z)~ 3 / 2 [q'^(z)]~ 1 / 2 is defined in terms 
of the luminosity distance di = c(l + z)Hq 1 J z dz'[q\ (z')]^ 1 , where q\(z) = 

^J{1 + Vt m z){l + z) 2 — Vt\{2z + z 2 ). Results presented here are for an Vt m = 
0.3, = 0.7 cosmology with h = 0.65. Figure 1 shows d L and d cosmo for this 
cosmology. 

The trigger efficiency P tr ((f) ph ) = exp[-((f) phfi /(f) ph ) a P h ] [27,7], where cf) ph is a pho- 
ton flux. On the 1024 ms timescale, <p p h,o — 0.26 ph cm -2 s^ 1 and a p h = 5.3. 
In the 50-300 keV band, the typical photon energy is ~ 100 keV when the peak 



of the vF v spectrum lies in the BATSE energy range. The vF v peak flux thresh- 
old is therefore PiO — 1-6 x 10~ 7 ergs x0 p?iiO ~ 4 x 10~ 8 ergs cm -2 s _1 . Thus 
P tr [(p p (z)] — >■ H[<f>p(z) — (f> P fl], where H[u] is the Heaviside function such that 
H(u) = 1 for u > and H(u) = otherwise. 

The flux density 5(e) = (1 + z)dJ j 2 dL[e(l + z)]/dfl, where the directional spec- 
tral power at frequency v = m e c 2 e/h is dL(e)/dQ. Although beaming can affect 
observable properties of a given GRB and the average event rate and total number 
of events in the Galaxy it has no practical effect upon average power, because 
the product of the number of sources and energy or power per source remains 
unchanged. Only apparent isotropic energy releases and powers are quoted hence- 
forth. 

The vF v flux threshold flux (f) p for a GRB with S(e) oc e _1 is given by the 
relationship <p p = {^Kd\Y x L peak = (47rd|) _1 L / ln(e max /e min ). As we argue below, 
detectors are most sensitive to GRBs when the photon energy e° of the peak of vF v 
flux is within the detector frequency window. In this approximation, the bolometric 
power L Q exceeds L peak by the bandwidth correction factor hi(e max / e m i n ) , which 
might be as large as ~ 5-10. 



A Clean and Dirty Fireballs 

Consider a GRB detector that is most sensitive in a range centered about dimen- 
sionless photon energy e. The peak luminosity of a fireball is [1] 

f td 1 ^)' 5 » dirt y fireballs with e > e° 
L P eak ^ rj E o < (3) 
*J 1 (-xr) 4/3 , dean fireballs with e « e° 

where e p is the photon energy of the peak of the vF v spectrum during the prompt 
7-ray luminous phase of a GRB, Eq is the total energy released in the explosion, and 
i] is the efficiency to transform this energy into radiation during the deceleration 
timescale [6] 



[l + z)x d 1 + z , 3E , 1/3 _ E, 



cIq cIq 47rp lo n 2 

Equation (3) refers to a uniform circumburst medium with rest mass energy density 
p = nm p c 2 , where the proton density n = 100n 2 cm~ 3 , E 54 = _E /10 54 ergs, and 
Xd is the deceleration distance. The efficiency factor rj that relates L peak with E 
and spectral parameters is [1] 

0/0 oA ( 0.03 , adiabatic limit with q = 1.6 

v - 3 ^ ~ ^ ~ ( 5 ) 
2^(3/4 + 5- 1 ) [ 25 ^ radiative limit with ^ = 2.9, 



where 2 + 5 is the photon number index of the prompt GRB spectrum at e £ e°, 
and the quoted efficiencies are for 5 = 0.2. The term g defines the radiative regime 
according to the asymptotic behavior of the blast wave Lorentz factor T(x) oc 
(x/x d )~ 9 when x ^> x^- For a roughly adiabatic blast wave with g = 1.7, as 
inferred from fits to the BATSE data [7], rj ~ 5-10%. 

In the external synchrotron-shock model [30,31], e° = qpy 2 Yy{l + z), where q is 
a parameter that takes into account the magnetic field energy density and efficiency 
to deposit swept-up energy into electrons. Assuming q is constant with respect to 
changes in To, £"0, and po, then we can define the Lorentz factor To such that the 
observed peak vF v frequency is equal to e . Hence gp^f q = e defines F , and we 
find that 

p l/3+<5/2 p8/3+45 ^ for dirty ^ ^ < f „ 

(6) 

p 1//3 r 8 ^ 3 , for clean fireballs with T ^ T . 

As can be seen from equation (6), there are enormous selection biases against de- 
tecting dirty fireballs because the peak luminosity measured by the detector varies 
£ Tq 5 when r £ f . A steep reduction in the peak luminosity at photon energy 
e likewise appears to be the case for clean fireballs, but the shorter durations of the 
clean events (see eq. [4]) mean that detectors trigger on fluence rather than peak 
flux when T 3> f . Thus GRB detectors are not so biased against the detection 
of clean as dirty fireballs, yet such detectors still preferentially detect those GRBs 
that produce prompt radiation with a vF v peak in the waveband at which the de- 
tector has greatest effective area. To avoid fine-tuning of T , one realizes that there 
must be two heretofore "invisible" classes of explosive phenomena that remain to 
be discovered by detectors with appropriate design [1]. These considerations also 
resolve the apparent paradox that a beaming scenario should produce a wide range 
of €p, contrary to observations [2,7]. 

B Comoving Rate Density of Fireball Transients 

Various lines of evidence lead to the conclusion that GRBs are formed in gaseous 
star-forming regions and are events related to the collapse of massive stars to black 
holes, possibly through the intermediate formation of a rapidly spinning neutron 
star (most similar to the supranova model [28]). Calculations of GRB statistics are 
simplified when an assumption is made that the comoving FT rate density h com (z) 
is proportional to the star formation rate (SFR) history of the universe [29] as 
traced, for example, by Hubble Deep Field and SCUBA data. 

As the universe ages, we can expect star formation to be triggered by galactic 
mergers and interactions and therefore to be proportional to the density of galaxies. 
In this simple way of looking at things, the star formation rate is thus proportional 
to the proper galaxy density n proper oc (1 + z) 3 n com near the present epoch. At 
earlier times, the galaxy formation rate would tend to decline due to the time it 
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FIGURE 2. Parameterization for comoving rate density of bursting sources for different values 
of parameters a and b, and with £ set equal to unity. The curve with a — 0.03 and b = 1 is 
similar to the star formation rate as traced by UV radiation in the Hubble Deep Field [32]. 



takes to form the first structures, which depends on the primordial spectrum of 
density perturbations and the physics of galaxy formation. Star formation would 
therefore be impeded at recent times and at early times. 

This suggests that we parameterize the comoving rate density of FT sources by 
the expression 



a) 



■ ( gfl 

ncom[Z) [(1 + z)~* + a(l + z)>] 



(7) 



which is normalized to £ at z — >■ 0. Figure 2 shows a range of parameterized models 
for the model SFR rate and, by assumption, the rate density history of explosions 
of a certain type. A model with a = 0.03 and 6 = 1 provides a reasonable facsimile 
of the SFR curve reported in Ref. [32]. 



C GRB Redshift and Size Distributions 

The left panel in Fig. 3 shows the observed redshift distribution of GRBs com- 
piled from Jochen Greiner's tabulations 2 . The right panel shows a model redshift 
distribution that is obtained by differentiating equation (2) with respect to z, giv- 
ing dN/dz = ATrcdl osmo (z)n com (z) / H . This model assumes that the rate density of 
GRBs follows the parameterization of equation (7) with a = 0.03 and b = 1, and 
that all events are detected. In the absence of data from an extremely sensitive 
GRB telescope, this model is not directly comparable with data. Nevertheless, it 
can be seen that both distributions peak near z ~ 1-2 and display a tail extending 
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FIGURE 3. (left) Distribution of GRBs with measured redshifts. Solid histogram represents 
all GRBs, and dotted histogram represents GRBs with Ej > 10 53 ergs, (right) Calculated dis- 
tribution of redshifts assuming that the comoving event rate density is proportional to the star 
formation rate history of the universe as parameterized by equation (7) with £ = 1CP 90 cm~ 3 
s _1 . Inset shows the low-z distribution. 

to high redshifts. The observations suggest an excess of GRBs with low redshifts in 
comparison with the model, and this feature would be noted even if GRB 980425, 
associated with the SBb host galaxy of SN 1998bw at z = 0.0085, were excluded. 
The redshift distribution of GRBs with z ^ 0.5 is posed as a crucial discriminant 
between models [3,7,25], because a large population of low luminosity, low-z GRBs 
(i.e., dN/dz ~ const for £ z <, 1) requires a source rate density ~ 2-3 orders of 
magnitude greater than would be the case if there were no low luminosity or low 
kinetic-energy GRBs. Although the FT rate would be much larger if a population 
of weak GRBs were discovered, the additional energy involved would be only a 
small fraction of the total energy released by GRBs with £"54 £ 10~ 2 [3]. 

To understand this, we plot size distributions for different peak luminosities L pea k 
in Fig. 4a using equations (2) and (7) with a = 0.03 and b = 1, using our chosen 
cosmology. These size distributions give the detection rate of GRBs that exceed 
peak vF v fluxes <p p , and are comparable to peak count-rate distributions reported 
by GRB detectors such as BATSE that trigger over a relatively narrow bandwidth 
when the count rate exceeds several a over background on some specified timescale. 
The bolometric peak luminosity L = BL pea k [23] for a spectrum uF u oc const, 
where the bandwidth correction factor B = \n(e max / e m i n ) ^ 1. Strong selection 
biases against detection come into play for sources that are detected at peak vF v 
frequencies e p which are much different than the detector's sensitive frequency e 
(Section IIIA). 

The size distributions approach the Euclidean behavior N(> 4> p ) oc (p p 3 ^ 2 at the 
bright end, and flatten to a constant at the dim end. The turnover flux in the size 
distribution occurs near the value of 4> p that corresponds to the peak vF u flux that 
would be detected from a source with brightness L peak at z = 1. Fig. 4b verifies 
this. 
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FIGURE 4. (a) Size distribution of GRBs as a function of vF v peak flux <j> p , for different peak 
luminosities L pea k and with S = 1CP 90 cm" 3 s" 1 . (b) Peak flux <p p as a function of redshift z for 
different values of L pea k- 



The redshifts of GRBs made possible by Beppo-SAX and the IPN (Interplan- 
etary Network) and GCN (GRB Coordinates Network) display a range of peak 
luminosities of GRBs spanning at least 3 orders of magnitude. The redshifts of 
these sources are distributed about 0.5 £ z p £ 2, though with unclearly defined 
wings at both low and high redshifts (Fig. 3a). To arrange a distribution of mixed 
peak-flux values in nearly the same redshift range requires a distribution dN/ dL pea k 
of L pea k values that vary £S L~^ ak with s R3 3/2. The pronounced change in slope 
of the peak count rate size distribution measured with BATSE near count rates of 
pa 8 ph cm" 2 s" 1 [33] or <p p R3 1CT 6 ergs cm" 2 s" 1 implies from Fig. 4a that the L peak 
distribution is flatter than —3/2 at L peak ^ 10 52 ergs s _1 and steeper than —3/2 at 
Lpeak ~ 10 53 ergs s -1 . The distribution cannot however be much flatter than —3/2 
at the faint end unless we dismiss the association of GRB 980425 with SN 1998bw, 
or at least relegate it to a separate class. But this would then compromise the 
arguments that equate reddened optical excesses observed in the late-time optical 
afterglows of GRB 970228 and GRB 980326 with SNe emissions [34-36]. 

Thus we assume that the peak luminosity distribution of GRBs follows a 
dN/ dL pea k oc L elk behavior over a wide range of L pea k- To the extent that L pea k 
is related to Eq through an efficiency factor rj, we can also parameterize the GRB 
E distribution by dN/dE oc Eq 3 ^ 2 . This behavior was originally derived from the 
numerical model [7] but lacks a fundamental understanding at present, though it 
clearly reflects an underlying convex function describing the distribution of GRB 
events in terms of apparent energy release. 



D Fireball Transient Rate in an L* Galaxy 



From Figs. 1-3 we see that the integrand in the calculation of the peak flux 
distribution in equation (2) is dominated by events near z p ~ 1-2, at least for 
sufficiently sensitive GRB detectors that reach threshold fluxes (j) Pi thr *C 1CT 6 ergs 
cm -2 s -1 . From equation (2), 

AlXC r z max{<frp) AlIC 

N (> <t>p)--FT dz rf cosmo(^) n com (z) ~ — -Az- d cosmo (z p ) ■ n com (z p ) , 

n Jo n 



which applies when the sensitivity threshold of the telescope is sufficiently good 
to detect sources at z ^> z p . This appears to be the case for explosions yielding 
7-ray energies E 1 ^ 10 53 ergs that represent about 1/2 of the GRBs with measured 
redshifts (see Fig. 3a). The mean redshift of this sample is ~ 1.5, and Az = 1-2. 
This allows us to derive an absolute normalization on n com (z p ) and therefore the 
present day comoving rate density of FTs that produce events with E y £ 10 53 ergs. 

From Fig. 1, we see that d cosmo (z p ) = 0.5 x 10 28 cm. Writing the observed burst 
rate in terms of the frequency ub of observed events per day, we find from equation 
(8) that 

n com {z p ) ~ 3.0 x 10- 90 ^ cm" 3 s' 1 . (9) 

The FT rate in an L* galaxy depends on its effective comoving volume. The 
density of L* galaxies in the local universe is ~ 3 x 10 6 Gpc -3 [3], so that an 
effective volume of an L* galaxy is Vl* ~ 10 76 V76 cm 3 , with V76 ~ 1. This number 
has a large uncertainty and could depend upon redshift. Thus the rate at which 
an L* galaxy produces FTs with E y £ 10 53 ergs is 

N*(E, > 10 53 ergs; z) = h com V L , = 3 x 10" 14 V -^L^ ~ 0.9^^ GEM , 

(10) 

where GEM is Galactic events per million years [37] . The factor 

T{z) = n com (z) (n) 

°^com\^" P ) 

adjusts the comoving star formation rate averaged over the redshift range extending 
to R3 2 x z p to the SFR rate occurring at a particular epoch z. To determine the 
event rate at the present epoch in the Milky Way, supposing it to be a typical L* 
galaxy, 

F ° - ^ - & 3m+t > + °& V( * Vi, a = o.o3 - °- 126 - (12) 



Note that the peak redshift z p of function (7) is z p = — 1 + (3/a6) 1 ^ 3+6 \ 

Two selection biases must be added into the accounting of equation (10) to 
complete the problem. One is to account for the dirty and clean fireballs that do 
not trigger GRB detectors. The dirty fireball contribution is probably the more 
numerous unseen population, and the GRB statistical analysis in Ref. [7] shows that 
dN/dT oc Tq 25 , for T <, 260, with dN/dT falling off more steeply at larger values 
of To- 3 A simple integration gives the detection fraction of detectable transients 
compared to all FTs with Tq > 2. If GRB detectors lose sensitivity to GRBs with 
T < 100 and 200, then the detection fraction is 0.5 and 0.18, respectively. This 
implies a clean and dirty fireball bias factor gr ~ 2-5. 

The second is to account for the environmental diversity that affects detection. 
Although blast wave physics shows that the emission properties scale to first order 
as ~ PoTq [1], the received flux gets washed out by light travel-time effects at low 
external medium densities, and by photoelectric absorption, Compton scattering 
and pair production attenuation [38,39] at high external medium densities. The 
low density bias seems unlikely given evidence that GRBs form in galaxies that 
are undergoing active star formation, but the high density bias is quite probable 
given recent results from Beppo-SAX [40] and Chandra [41] that seem to require 
a highly metal-enriched medium from SN events that precede by some months to 
years a GRB in the cases of GRB 990705 and GRB 991216, respectively. This is 
rather difficult to quantify, but an environmental biasing factor might be in the 
range 1 £ g po £ 10, implying a total biasing factor g = Qv a Q P0 ~ 2-30. 

In view of these considerations, we find that events with > 10 53 ergs occur in 
a local L* galaxy within z £ 0.1 at the rate 

N^E, > 10 53 ergs) (GEM) = ^s^og = .ty . (13) 

Az 

Our best guess is that ip ~ 1. Equation (13) sets the normalization for the event 
rate of FTs in an L* galaxy such as the Milky Way. 



E Power and Number of FTs in the Milky Way 

GRBs detected with energy E 1 are produced by explosions with total energy 
E = Ej/rj ~ 10-100 E 7 , using equation (5). We have argued in section IIIC that 
the E'o-distribution of GRBs can be approximated by 

N*(E ) = ^(|^)~ 3/2 (14) 

where the exponent —3/2 is chosen in view of the observations which indicate a 
wide range of values of E y at z ^ 1-2. We now adopt a notation where £ refers to 

3 ) The density of the uniform external medium was set equal to 100 cm~ 3 , but the spectral model 
[1] is degenerate in the quantity uoTq. 



energy in units of 10 54 ergs. Normalizing equation (14) to the L* rate of GRBs and 
FTs with £ < £ < E\ implies that 



K = K(r ] ,£ ,£ 1 ) = )l . (15) 

Z[C — C-i \ 

It now only remains to assign £\ and relate £q to £ 1 for a given efficiency rj, which 
in turn relates N*(Eq) to iV*(£ 7 > 0.1) from equation (13). 

The efficiency rj to convert the injected kinetic energy E into 7-ray energy E 7 , 
given by equation (5), is assumed to be constant from event to event. We consider 
three efficiencies, namely rj = 100%, 10%, and 1%, which imply a maximum value 
of £\ given the beaming factor. The maximum energy available from the collapse 
of a neutron star, even a rapidly rotating 3 M neutron star, is ~ 10 54 ergs. 
A beaming factor SQ/An £ 1% seems ruled out by the data for the onset of the 
beaming break in some GRBs such as GRB 990510 and GRB 990123, as well as 
by the interpretation of fluorescence Fe Lya emission [41]. Thus we assume that 
the maximum apparent energy release is £\ = 100, corresponding to a 1% beaming 
factor, and £2 = 10, corresponding to a 10% beaming factor. 

Our grid of model runs therefore contains Runs 100a, 100b, 10a, 10b, and lb, 
where the number in the label refers to the percentage value of 77, and the letters 
a and b refer to £\ = 10 and 100, respectively. £q takes values of 0.1, 1, and 10 for 
rj = 100, 10, and 1%, respectively. There is no Run la because £q = 10 = £\. 

From K, we obtain the total power of FTs into the Galaxy from the expression 

r£i 

= L*(> £ ) = KE 54 / d£ ■ £- 1/2 ~ 2 x lO 54 /^ 72 ergs s" 1 (16) 
for £ <C £1, and we obtain the total number of FT events from 

iV*(> £) = tgaiK I" 1 d£' £'- z ' 2 ~ 2t gal K£- 1 ' 2 (17) 
J s 

for £ <^ £1, where the age of the Galaxy is given by t ga i = 10 10 ti yr. 

Table 1 gives the results of K, the FT power L*, and the total number of FT 
events TV* during the life of the Milky Way for this grid of models. The value of 
iV* is sensitive to £ which ranges from 10~ 6 ^ £ £j 10~ 3 , corresponding to events 
with 10 48 £ .Eo(ergs) £ 10 51 . The lowest value of £ seems required if GRB 980425 
is believed to be associated with SN 1998bw, implying a value £ £ 10~ 3 . 

As can be seen from Table 1, FTs inject a time- and volume-averaged kinetic 
power between ~ 10 39 -10 41 ergs s _1 into the ISM, depending sensitively on the 
assumed efficiency for converting the directed kinetic energy of the outflow into soft 
7 radiation. Comparing with equation (1), we see that FTs inject adequate power 
to produce the hadronic cosmic rays if the CR hadrons are accelerated by GRB 
blast waves with high efficiency, yet radiate gamma-rays with modest (rj ~1-10%) 
efficiency. A value of i] ~ 5% is implied by calculations from the external shock 



TABLE 1. FT Power and Number of FT Events in an L* Galaxy such as the Milky Way 
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model [7]. A value of r\ ~ 1%, as seems appropriate for an internal shock model, 
would require more power from GRB progenitors. However, no good statistical 
treatment has been performed within a wind scenario that can be used to assess 
detector biases, and these could change the value of ip. For example, the beaming 
paradox resolved in Section IIIA in the external shock model remains unsolved to 
date by an internal shock model. But in either case, the progenitor stars of GRBs 
are seen to inject £ 10 40 ergs s" 1 into the ISM of an L* galaxy like the Milky Way. 

IV COSMIC RAYS, SUPERNOVAE, AND FIREBALL 

TRANSIENTS 

The most crucial test for a theory of cosmic-ray origin is to demonstrate that 
the putative source class provides the available power. The much rarer FTs and 
GRBs can provide this power because they are also much more energetic. An 
important point is that SNe are thought never to liberate more than a few xlO 51 
ergs of ejecta kinetic energy per explosion, with the bulk of the power that is 
potentially available in core-collapse SNe carried away by neutrinos. The massive 
stellar collapse events that make FTs evidently liberate as much as 10 54 ergs per 
explosion in the supranova model; perhaps an order-of-magnitude more energy is 
possible in a collapsar model [42] or in the collapse of a massive Fe core to a 
black hole [3]. More likely, there is no strict uniformity of the properties of stellar 
progenitors that collapse to form black holes. 

A Ultra-High Energy Cosmic Rays 

The energy density of UHECRs with particle energy £ 10 20 eV (see Fig. 9 in 
Ref. [3]) is 

u v ~ few x 1(T 21 ergs cm" 3 ~ 77- Vy (18) 

where r P7 = A P7 /c is the timescale to degrade a particle's energy by a factor of ~ 2 
through photopion production of UHECRs with the cosmic microwave background 



radiation. The term L*/V* is the average power injected into an L* glaxy divided 
by the effective volume of an L* galaxy. For 10 20 eV protons, the mean energy loss 
length is ~ 140 Mpc [43]. Letting = 10 40 L 40 ergs s -1 , and recalling that the 
effective L* galaxy volume V ga i = 10 76 Viq cm 3 , we obtain 

^r P7 ~ 1.4 x 1(T 20 ^° ergs cm" 3 . (19) 

If FTs and GRBs power the UHECRs, some modest but nonnegligible factor £ 0.01 
of the available FT and GRB power must be processed into UHECRs to account 
for the observed energetics of metagalactic UHECRs. Waxman [44] and Vietri [45] 
first pointed out this coincidence. 



B Cosmic Rays between the Knee and the Ankle 

The energy density of cosmic rays above the knee at = 2-4 PeV is 

Mknec ~ 10" 16 ergs cm" 3 . (20) 

If GRBs and FTs power the PeV - EeV CRs, then CRs above the knee of the 
spectrum diffuse from the disk into the halo and then escape into intergalactic 
space. 

Rather than attempt to derive the halo escape timescale th from a model for the 
Galactic halo, we instead invert the relation 

^r H ~ 10~ 16 ergs cm" 3 . (21) 

to solve for t h . Let the halo volume V H = nr 2 £ and the escape timescale r H = 
kni/c, where r = 15 kpc is the galaxy radius and £ is the effective halo height. 
Setting -u kne e c± L*kH/(nr 2 c) = 5 x 10 -17 L^kn ergs cm -3 , we obtain 

L m k H ~ 10 . (22) 

This relation implies interesting properties concerning the halo diffusivity. If L 40 ~ 
1 for the power into the PeV- EeV component, then the escape time of CRs near 
the knee of the spectrum is only an order-of-magnitude greater than the transit 
timescale, implying a weak halo magnetic field. 

The high-mass progenitor stars of FTs could therefore power the PeV- EeV CRs 
if £ 10% of the energy of FTs is transformed into particles with these energies. 
The bend at the knee of the CR spectrum is due, in this interpretation, to galactic 
propagation effects. 



C GeV - PeV Cosmic Rays 



The galactic cosmic rays between 10 9 -10 15 eV require, if distributed uniformly 
throughout the disk of the Galaxy, a power ~ 5 x 10 40 ergs s _1 (Section II). For 
efficiencies ~ 1-10% to transform the explosion energy into observed 7 rays, Table 
1 shows that FTs and GRBs in our Galaxy could satisfy this power requirement, 
though the efficiency to transform the directed kinetic energy into CR energy must 
be high. Stars that collapse to black holes could therefore be the sources of the 
hadronic cosmic rays. 

The argument that the collapse of GRB and FT progenitor stars to black holes 
accelerates both the UHECRs and the GeV-EeV hadronic CRs is primarily based 
on simplicity — this solution requires the fewest number of theoretical assumptions. 
A SN origin for the GeV- PeV Galactic CRs leaves open the question of cosmic rays 
with energies above the knee of the CR spectrum; FTs as the sole source class 
closes it. Compositional and theoretical studies near and above the knee of the CR 
spectrum will provide much of the necessary data to answer it. On the basis of 
available power from these explosions, it is therefore proposed that the high mass 
stars that collapse to form GRBs and FTs are the sources of the hadronic CRs. 

Other arguments for this cosmic-ray origin hypothesis are given in my paper 
written in memory of Jan van Paradijs [3]. 

V THE UNIDENTIFIED EGRET SOURCES 

Apparently independent of the GRB problem looms the question of the identity 
of the ~ 15 COS B, > 100 EGRET, and soon-to-be ~ 10 3 (?) GLAST unidentified 
sources. These objects comprises an ~ 3° half-angle scale height disk component, an 
~ 20° half-angle scale height mid-latitude population (which is wider and more ex- 
tended than a bulge population), and an isotropic component representing £ 10% 
of the unidentified EGRET sources (UES). There is also detected a 10° diameter 
crown of a dozen or so EGRET 7 ray sources centered at I — 0°, b — +5°, that is 
incidentally similar to a feature in COBE/DIRBE maps. 

Undoubtedly, some of the UES are pulsars, given our knowledge of Geminga. 
Grenier [46] presents compelling evidence that some of the disk sources are coinci- 
dent with X-ray plerions, which opens the possibility that several other low-latitude 
sources are members of this class. There may be separate populations of disk and 
mid-latitude sources [47], with the remaining isotropic fraction being some uniden- 
tified blazars of an extreme type, for example, COMPTEL MeV blazars or extreme 
blazars that put the vF v peak of the nonthermal synchrotron emission in the soft 
7-ray band [48]. Some of the isotropic component could consist of the synchrotron 
self-Compton emissions of dirty fireballs [51] caught once in the act. Here the better 
imaging of GLAST will be enormously helpful. Each EGRET unidentified source 
is, of course, potentially a TeV source. 



On the basis of the spatial, spectral, and variability data, I [52] concluded that the 
solution to the UES which entails the least number of unproven assumptions is that 
these sources are isolated black holes accreting from the ISM. Only one population 
of such sources is required; the different scale-heights of the low-latitude and mid- 
latitude components represent younger (~ 10 7 -10 8 yr old) black holes found in 
dense molecular cloud complexes, and older black holes with large scale heights 
that are located within several hundred pc of the Solar system, respectively. 

The spectral and variability data argue against a large number of pulsars or 
isolated neutron stars comprising the UES if their emissions are like Geminga. The 
UES have soft spectra like blazars (photon index a ~ 2.3 ± 0.4) rather than hard 
spectra like pulsars (a ~ 1.7, extending from 1.2-2.2). This argument cannot be 
pushed too far, though, as an underlying microquasar population has not been 
detected, and the spectral behavior of an isolated accreting black hole could be 
completely different from that of a blazar or pulsar. Bondi-Hoyle accretion from the 
ISM onto ~ 30-100 M Q black holes gives the correct order-of-magnitude luminosity 
for the ~ 10 35 ergs s _1 disk component and the ~ 10 33 ergs s _1 local component 
[52]. Ref. [53] points out that the luminosity relationship in advection-dominated 
models can alter the dependence on black-hole masses and Eddington rates, and 
advances this interpretation of the UES further. Punsly [54] and G. Romero, at 
the Alicante INTEGRAL workshop, have a different approach based on magnetized 
black holes. 

A model based on isolated black holes accreting from the ISM requires a few 
xl0 5 -10 6 black holes with masses between ~ 10 and 100 M & . Although this point 
requires more study, to satisfy compactness, variability and energy constraints, it is 
nearly certain that FTs with £ ^* 10 _1 entail black hole formation. Whether the 
FTs with Eq <C 10 _1 form black holes is not clear, but if they do, then as many as 
10 6 -10 7 black holes are scattered throughout L* galaxies in an irregular distribution 
reflecting their birth and later movements. 

Table 1 demonstrates that the requisite number of black holes can be formed. 
A population of £ 10 M black holes is formed in collapsar-type models [55,56]. 
The black holes produced in the supranova model [28] would have a few Solar 
masses unless delayed fallback of the ejecta onto the collapsing neutron star caused 
the black hole to become much more massive. Moreover, a few Solar mass black 
hole found in the vicinity of its own ejecta from an earlier SN event could grow 
quite massive if the ejecta has not yet been fully dispersed. In a region with 
effective neutral hydrogen columns exceeding ~ 10 25 cm~ 2 — and assuming that 
the subsequent GRB did not completely disrupt the dense ejecta — a nascent black 
hole could form an IR hotspot. Thus it may be possible even to have a population 
of moderate mass black holes in a supranova model. Whether this is consistent 
with IR source counts also requires further study. 

It hardly needs stated that AGILE, GLAST, and the new generation of TeV 
telescopes will make enormous progress in characterizing spectral, variability, and 
counterpart behavior of UES. Microlensing of background radiations to detect iso- 
lated black holes [49] looks, upon initial inspection, extraordinarily difficult with 



present technology. Binary black hole systems such as Cygnus X-l are well-known 
[50]; GRBs might announce their birth. 



VI THE LINEAGE OF STARS 

The number of SN events that have occurred in the Milky Way is ~ 10 10 yr/30 
yr i=s 3 x 10 8 and, as we have seen, FTs leave behind ~ 10 4 - 3 x 10 6 black holes 
ranging in mass from a few Solar masses to tens of Solar masses. From what range 
of stellar masses do FTs originate? 

TABLE 2. Rates of Different Nova Types in Supernova Units 





Nl 


N2 


N3/N4 


NF & 


Sbc-Sd 


0.21 


0.86 


0.14 


0.003 



a Computed in [3] . 

This question can be answered by considering explosive event rates in Type Sbc- 
Sd galaxies, using the tabulation of Ref. [57]. (The Milky Way is Type Sbc, though 
with a bar of uncertain size.) The Supernova Unit is the number of events per 10 10 
Solar blue luminosities Lb,q per century. If the Milky Way radiates 2xl0 10 Lb,q, 
then Table 2 implies a time-averaged rate for all SNe types of 1 per 60 yr. Note the 
monotonic trend in Table 2, with the exception of the Nls (SN la's). The Nl rate 
is below the trend in the rates of core collapse novae because white dwarfs depend 
in most cases on being formed in binary systems to be driven into detonation. The 
ratio of the rates of N3/4s (SN Ib/Ic's) compared to the rate of all core collapse 
events (N2, N3, N4, and NF) is ~ 0.14. The ratio of the rates of NFs (GRBs 
and FTs) to all core-collapse events is ~ 5.0 x 10~ 3 iV 6 for a 10 10 yr lifetime of the 
Galaxy, assuming that FTs make lO 6 ^ black holes during this period. 

The differential stellar IMF (initial mass function) £(JA) oc A4~ a , where Ai is 
the mass of the zero age main sequence (ZAMS) star in units of Solar mass, and 
o is the differential IMF index. For a Salpeter IMF, a s = 2.35, and a st = 2.8 for 
a steep IMF. If Mthr is the mass threshold for initiating core collapse events, then 
the fraction of stars with mass > M. is [58] 

- a hr dMaM)- - Uand5x10 N ^ (23) 

corresponding to the ratio of the rates of N3/N4s and NFs, respectively, to the 
rates of all core-collapse events. 

The relation M. = M. t hr^ 1 ^ 1 ^ a gives the lower limit to the ZAMS mass of stars 
that eventually collapse as different types of novae. The mass threshold for making 
SNe II/N2s is argued [58] to be in the range Mthr — 6-14, depending in general 
upon metallicity. Stars with ZAMS masses 

_ f A.3M thr , a s = 2.35 , 
M > N ^ = U.0M thr , a st = 2.8 ' ^ 



i.e., with masses £ 20-60 M , end their lives as N2s and N3s. According to this 
simple argument, stars with ZAMS masses 

M ^ j51M thr N 6 - 74 , a s = 2.35 (t>t .. 
M>NF ~ { 19M thr N^ , ast = 2 .8 ' (25) 

i.e., stars with masses £ 110-300 M , end their lives as FTs. Similar results are 
obtained from the numerical calculations [3,7]. 

Based on these considerations, Fig. 5 sketches the scenario from the birth of 
stars on the stellar IMF with different ZAMS masses to their final fate as novae of 
various types, including the detonation supernovae, the NS core-collapse SNe, and 
the BH core-collapse hypernovae, namely GRBs and FTs. Observations and stellar 
structure calculations will test this picture, it being recognized that the actual 
situation is far more complex. See Ref. [3] for more discussion of this classification 
scheme and implications for different types of galaxies. 



VII COSMIC RAYS FROM GAMMA-RAY BURSTS 

The statistical estimates of available energy and power precede all subsequent 
considerations relating to a cosmic-ray source model. This contribution shows that 




FIGURE 5. Sketch indicating the relationship between the stellar IMF £(M), the ZAMS mass 
of the progenitor stars, and the explosion types that result from stars with different ZAMS masses. 
The shaded area represents the low ZAMS mass stars that evolve into white dwarfs on time scales 
exceeding the age of the Galaxy. 



the available power from the high mass stars that collapse to form GRBs and FTs 
is sufficient to power cosmic rays. 

Space limitations do not permit me to address the cosmic-ray acceleration and 
adiabatic loss problems here. Because they were such a strong focus of criticism 
of this CR origin hypothesis at the Heidelberg 7 2000 workshop, however, it seems 
worthwhile to outline a proposed solution. 

Future writings [12] will demonstrate that gyroresonant stochastic particle ac- 
celeration from turbulence in the relativistic blast waves of GRBs and FTs can 
accelerate protons and ions to ultra-high energies; that stochastic acceleration of 
particles by the evolving turbulence spectrum preferentially accelerates particles 
with a number index somewhat steeper than ~ —2; that £ 10 19 eV protons and 
ions can leak out of the blast wave during the prompt and afterglow phases to 
form the UHECRs, their Larmor radii being so large that they avoid adiabatic 
losses; that the bulk of the accelerated particles remain trapped in the blast wave 
throughout the prompt, afterglow, and nonrelativistic Sedov phase during which 
shock Fermi acceleration becomes more efficient than stochastic acceleration; that 
the late phases of a FT/GRB remnant represent in many respects a SNR except 
that higher energy particles are available from the acceleration that takes place as 
the relativistic blast wave decelerates; and that cosmic rays finally leave the rem- 
nant as it dissipates in the ISM, so that the adiabatic loss problem is solved (as 
in the standard model) by shock acceleration taking place during the expansion of 
the SNR until the energy density of the relativistic particle fluid is small compared 
to the magnetic field energy density of the ISM. 

Observations may well refute the cosmic-ray origin hypothesis put forward here. 
GLAST observations of a distinct ir° feature in the vicinity of young SNRs could 
provide the 7-ray evidence that has been sought for decades. TeV telescopes could 
detect emission from SNRs that is more likely to have a hadronic than leptonic 
origin. On the other hand, GLAST or TeV observatories could detect hadronic 
emission at sites that arguably witnessed an earlier GRB, such as the Cygnus 
region, the Sco-Cen complex, or at sites in the direction of the ~ 10 18 eV cosmic 
ray anisotropies discovered by the AGASA [59] and SUGAR arrays [60]. 

In summary, the available evidence is not yet adequate either to verify or to 
reject a SN or GRB/FT origin of the hadronic cosmic rays. 
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